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respect to the Bures measure and show that it is equal to that of a A''^ — 1 dimensional hyper- 
hemisphere of radius 1/2. For N — 2 we obtain the volume of the Uhlmann hemisphere, ^S'^ C R"*. 
We find also the area of the boundary of the set Mn and obtain analogous results for the smaller 
set of all real density matrices. An explicit formula for the Bures-Hall normalization constants is 
derived for an arbitrary N. 
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I. INTRODUCTION 

Modern applications of quantum mechanics renewed interest in the properties of the set of density matrices of finite 
size. Its geometry depends on the metric used 0,12 • However, the simplest possible choices exhibit some drawbacks: 
the trace metric, defined by the trace distance 

D^,{P,<^) = Hp- 'Till :=Tr|p-a| =TrV(p-^)2, (1.1) 

is monotone but not Riemannian, while the Hilbert-Schmidt metric, induced by the Hilbert-Schmidt scalar 
product, {A\B) = TiA^B and related to the H-S distance, 

Dus{p,<y) = \\p-<7\\2--= v/Tr[(p-a)2] , (1.2) 

is Riemannian but not monotone . This fact provides an additional motivation to investigate other possible metrics 
and distances in the space of mixed quantum states. 

Nowadays it is widely accepted that the Bures metric 0,01, is distinguished by its rather special properties: it is 
a Riemannian, monotone metric, which is Fisher-adjusted [ij: in the subspace of diagonal matrices it induces the 
statistical distance 0- Moreover, it is also Fubini-Study-adjusted, since it agrees with this metric at the space of 
pure states . The Bures distance between any two mixed states is a function of their fidelity 0, 0| and various 
properties of this distance are a subject of considerable interest (see |i|l3,lil|ll|ll|ll|lll3). 

A given distance between density matrices induces a certain measure in this space. The probability measure in 
the space of mixed quantum states induced by the Bures metric was defined by Hall [3| • He demonstrated that the 
eigenvectors of a random state of size N generated according to this measure are distributed according to the Haar 
measure on U{N), derived an explicit expression for the joint probability density for eigenvalues of a random state 
generated according to this measure and found for = 2 the normalization constant Cn for the joint probability 
density of the eigenvalues. Further Bures-Hall normalization constants C3 and C4 were calculated by Slater but 
the problem of finding the constants Cpj for arbitrary A^ remained unsolved. 

The aim of this work is to compute the volume of the set at of all density matrices of size N and the hyperarea 
of its boundary according to the Bures measure. This work is complementary to our parallel paper |19| |. in which the 
volume was computed according to the Hilbert-Schmidt measure, but the notation and the structure of both papers 
are not exactly the same. In this work we treat in parallel two different cases: of complex density matrices, and its 
particular subset, of real density matrices. Both cases are labeled by the 'repulsion exponent' (3 = 2 for the complex 
case and /3 = 1 for the real case. 

The paper is organized as follows. In section II the Bures distance is defined and its remarkable properties are 
reviewed. In section HI we derive the Bures measure and express the volume and the surface of the space of mixed 
states as a function of the Hall normalization constants. Key results of the paper are presented in section IV, in which 
we compute the normalization constants Cn for arbitrary N, see (|4.11ll . and calculate explicitly the Bures volume 
of A^AT, see Eq. (|4.12|) . In section V we perform a similar analysis of the space of real density matrices. A concise 
presentation of the results obtained is provided by Eq. H5.15|) , from which one may read out the volume of the sets of 
complex and real density matrices, the surface of their boundaries and their 'edges'. 
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II. BURES METRIC 



We shall start the discussion considering the set of pure states belonging to an TV-dimensional Hilbert space 
TCn- This set has the structure of a complex projective space, CP^~^, and there exists a distinguished, unitarily 
invariant measure, which corresponds to the geodesic distance on this manifold. It is called Fubini-Study distance 

DpsilfA), I^b)) — arccos(2K — 1) = 2arccos(v^), (2.1) 

where k = \ {ipA\ipB)\'^ denotes the transition probability. This quantity was generalized for an arbitrary pair of mixed 
quantum states by Uhlmann 



F{pa,Pb) = [Tt{^pb ^/pI)'^'J = ll^/pI^/p^ll^ (2.2) 

and its remarkable proprieties were studied by Jozsa 9] who called it fidelity |53j . It is a symmetric, non-negative, 
continuous, concave function of both states and is unitarily invariant. For any pair of pure states the fidelity reduces 
to their overlap, F = TrpAPB = |(V'a|V'-b)P = Hence a function of fidelity, called Bures length ^ or angle p^ . 

Da{pa, Pb) := arccos^F(/9A,PB) ]^&yccos{2F [pA, Pb) - l) , (2.3) 

coincides for any pair of pure states with their Fubini-Study distance (|2.1|l . 

Fidelity may also be defined as the maximal overlap between the purifications of two mixed states, |($A|'i'B)P, 
where |<i>i) S 'Hn^ — T^n T^n and both mixed states are obtained by partial tracing over the auxiliary subsystem 
E, Pi — Tr£:|$i)(4>i|. Fidelity allows one to characterize the "closeness" of a pair of mixed states: it is equal to unity 
only if both states do coincide. Another function of fidelity. 
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Db{pa,Pb) = {2-2^F{pA.pB)f'^ = ^2-2Tr^VP^PBVP^ (2.4) 

satisfies all axioms of a distance and is called Bures distance 0. Hence Bures metric is Fubini-Study adjusted, since 
for pure states it induces the same geometry. Several relevant papers on this subject are due to Uhlmann |^, I3, lioll . 
who showed that the topological metric of Bures is Riemannian. Hiibner provided a detailed analysis of the N = 2 \V\\ 
and = 3 case while Dittmann analyzed differential geometric properties of the Bures metric 12^. Relations 
of the Bures metric to information theory can be seen in papers by Fuchs and Caves j23], Vedral and Plenio 
while the closely related fidelity is often used in the comprehensive book by Nielsen and Chang [22| ■ It was shown by 
Braunstein and Caves 7] that for neighbouring density matrices the Bures distance is proportional to the statistical 
distance introduced by Woottcrs 25] in the context of measurements which optimally resolve neighbouring quantum 
states. 

For diagonal matrices, pA = diag(a) and pb = diag(6), the expression for the Bures distance simplifies and becomes 
a function of the Bhattacharyya coefficient [S^, B{a, b), 

N 

DUpA,PB) = 2-2j2\/aA = 2-2B{d,b) . (2.5) 

i=l 

This is the Euclidean distance between two points on a unit circle, the corresponding vectors form the angle 9 such 
that cos9 = a ■ b. The geodesic distance along the circle i? = 2arccos(l — Dg/2) — 2Da, is called in this case 
the statistical distance, since it provides a measure of distinguishability between both probability distributions by 
statistical sampling. Another function is called the Hellinger distance 



Dnia, b) 



\ 



N 



J2iV^^-Vhy^\/2-2Bid,b). (2.6) 



It is clear that for any two diagonal density matrices their fidelity reduces to the squared Bhattacharyya coefficient, 
F{p, a) = B'^{a, b), while the Bures distance between both states coincides with the Helhnger distance between both 
vectors of eigenvalues, Db{p, <j) = -Off (a, b). 
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a) Hilbert-Schmidt 



b) Bures 




FIG. 1: Geometry of the set of the mixed states for N — 2: a) - the 3-baU embedded in R induced by the Hilbert-Schmidt 
distance, b) - hemisphere of S'' embedded in R^, (shown as a cross-section of taken at y = 0), induced by the Bures distance. 
The maximally mixed state p« is located at the center of the ball (a), and at the hyper-pole of S"^ (b). The set of pure states 
is represented by the Bloch sphere (a) and by the Uhlmann equator £ (b). The cross-section of M2 with the plane 
z = produces a full circle in the equatorial plane (a), or a half of a hyper-meridian joining the states pw and pE (b). 



Let us write the Bures distance for infinitesimally closed diagonal states defined by the vectors a and a' — a + da, 

dsl = {DB{pA.p'A)f=-i-'^jZ^.\h + — ■ (2.7) 

, V li 

Since X^i^i = 1 X^i^i '^^i ~ ^' computing the Bures line element ds^ we need to expand the square root 
up to the second order, and obtain 

d4 = jE^=Tr(dV^r (2.8) 

- precisely the metric on a sphere, or also, the statistical metric on the simplex, also called Fisher-Rao metric psj . 
Hence the Bures metric is Fisher-adjusted - restricted to the simplex of diagonal matrices it agrees with this metric. 
For N = 2 the simplex of diagonal matrices is one-dimensional, and from the point of view of Bures metric it forms 
a semicircle. 

It is instructive to compare the geometry induced by the Bures (|2.4ll and the Hilbert-Schmidt metric (|1.2|l . Any 
N — 2 density matrix p may be expressed in the Bloch representation 

p=i(l + f.a) (2.9) 

where a denotes a vector of three Pauli matrices, {ax,(Jy,(Jz}, while the Bloch vector t = {x,y,z} belongs to M'^. 
With respect to the HS metric the set of mixed states A42 forms the Bloch ball with the set of all pure states 
at its boundary, called Bloch sphere. The H-S distance between any mixed states of a qubit is proportional to the 
Euclidean distance inside the ball, so the H~S metric induces the flat geometry in A^2- The radius of the Bloch ball 
depends on the normalization of the Bloch vector, and is equal to unity while using the definition (|2.9|) . 
Working in a Cartesian coordinate system we may rewrite ()2.9|l as 



1/ 1 + z X — iy 
2 \ X + iy 1 — z 



x^ + y^ + z^<l (2.10) 



Setting Pa — p, Pb = P + dp in the expression for the Bures distance H2.4|l and expanding to second order we obtain 
one-fourth times the round metric on the 3-sphere, 

,2 1 /^j 2 , , 2 , , 2 , ixdx + ydy + zdz)'^\ 

dss = 7 +dy +dz +— —\ . 2.11) 
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Uhlmann noticed this fact in [lOj and suggested embedding of the three dimensional space of the N = 2 mixed states 
in M.^. The transformation 



p{x, y,z) - i^x, z, u := - - y2 _ z'^j (2.12) 

blows up the 3-ball into the Uhlmann hyper-hemisphere of S'^ of radius i?2 = 1/2. The maximally mixed state p*, 
mapped into a hyper-pole, is equally distant from all pure states located at the hyper-equator - see Fig. ^ The 
hyper-equator is equivalent to an ordinary two dimensional sphere and hence isometric to CP^. The comparison 
between both metrics is very transparent for real density matrices; the space of all rebits, which in the HS-metric has 
the structure of the 2-ball, in the Bures metric has the from of a hemisphere of S^. 

Since the radius Rb of the Uhlmann hemisphere is equal to 1/2 we obtain immediately the Bures volume V and 
the Bures area S" for = 2 

^ = voIb{M2) = ^SsR^B = ^ , Si^^ = vols (9X2) = S2R% = n, (2.13) 

where Sk denotes the volume of the unit sphere S''. This result was obtained by Caves ^2^, who also expressed the 
Bures volume of the set of mixed states in form of an integral. In the further sections of our work we are going to 
compute such integrals to derive an explicit result for the Bures volume for an arbitrary matrix size N. 

The last property of the Bures metric, we are going to review, concerns monotonicity. In the classical case of 
commuting diagonal density matrices, one studies the space of probability vectors a = {ai, ...ajv}, where > and 
^iLi = 1- This space is isomorphic with an — 1 dimensional simplex Ajv-i- A discrete dynamics in this space, 
a' — Ta, is given by an arbitrary stochastic matrix T of size N. This matrix contains non-negative entries, > 0, 

and the stochasticity condition, X^jLi = 1 for i = 1,...,N assures that the output vector a' is normalized. A 
distance in the space AAr_i of A^-dimensional probability vectors is called monotone, if it does not increase under the 
action of any stochastic matrix T, 

I)mon(a,&) > A„o„(r5,r6). (2.14) 

If a geodesic distance in the probability simplex is monotone, the corresponding Riemannian metric is called monotone. 
Chentsov showed '30'| that every monotone Riemannian distance in this space is a function of the Fisher statistical 
distance |2^, or the Bhattacharyya coefficient B defined in H2.5|l . 

The quantum case is more complicated. Any discrete dynamics may be described by a quantum operation, <I> : 
Mn ^ -Mn, which is trace preserving, Tr$(p) = Trp = 1, and completely positive. Positivity means that ^(p) > 0, 
i.e. $ sends positive operators into positive operators. From a physical point of view it is necessary to assure that an 
operation is well defined on the system coupled to an ancilla - an auxiliary system which describes the environment. 
Hence one requires that $ is conipletely positive, which means that any extension of the map by an identity operator, 
$(8) I, is positive (sj- Completely positive trace preserving maps play the role of stochastic matrices T in the quantum 
case, and are thus called stochastic maps. In full analogy to (|2.14|) any distance D in the space of mixed quantum 
states A^^^) is called monotone if it does not grow under the action of a stochastic map, 

Anon(p,CT)> i?mon($(/5),*(^T)) • (2.15) 

If a monotone distance is geodesic the corresponding metric on A4n is called monotone. 

In contrast to the classical case, there exist infinitely many monotone Riemannian metrics on the space of mixed 
quantum states. However, as shown by Morozova and Chentsov and Petz and Sudar 0|, they may be characterized 
by a symmetric function c{x,y). In general, the length of a vector (matrix) Bij may be position dependent; at a 
diagonal point A — diag(ai, 02, ajv) the squared length of B is given by P| 

N N 

II^IIa = E^ + 2 E c(a„afe)|B,fer (2.16) 

j = l ^ j<k=l 

For a diagonal matrix B the second term vanishes and the first term gives back the unique statistical metric on the 
simplex. The off~diagonal elements, which enter the second term in H2.16|l . allow for more freedom in the choice of 
a monotone Riemannian metric. The Morozova-Chentsov function satisfies c{sx,sy) — s~^c{x,y). The symmetry of 
this two-argument function c enables to express it by a function of a single parameter, c{x,y) = f{x/y)/y, such that 
f{l/t) = f{t)/t. 

A function / is called operator monotone, if for any matrices G, H the relation < G < H implies < /(G) < f{H). 
Here G < H denotes that their difference is semi positive, H — G > 0. Taking for G and H real numbers we see 




FIG. 2: Morozova-Chentsov functions f(t): minimal function for maximal metric (dotted line), Kubo-Mori metric (dashed 
line), and maximal function for the minimal, Bures metric (upper line). These metrics are Fisher-adjusted, /(I) — 1, while 
Bures metric is also Fubini-Study adjusted, /(O) = 1/2. 



that operator monotonicity implies that / maps R4. M+ and is concave. The result of Morozova and Chentsov 
[3^ , completed and extended by Petz and Sudar 0, Is^] , states that every monotone metric may be written in the 
form (|2.16() up a proportionality constant. Conversely, this equation determines a monotone Riemannian metric, if 
the function f{t) — l/c(t, 1) is: a) operator monotone, and b) self-transposed, f{l/t) = f{t)/t. Such functions / are 
called Morozova-Chentsov functions. Comparison of different metrics is meaningful under a certain normalization. 
Usually one requires: c) /(I) = 1, since for diagonal matrices such a metric agrees with the statistical (Fisher) metric, 
and is called Fisher-adjusted. There exists infinitely many functions / satisfying all three conditions a)-c), and 
leading to Fisher-adjusted metrics Let us mention three examples here: 

/minW^T^, /KM(t) = ^, /max(<) = ^ (2.17) 

t + 1 int 2 

containing both extremal cases and the intermediate Kubo-Mori (Bogolioubov) metric used in quantum statistical 
mechanics (34j and distinguished by the fact that both affine connections with respect to this metric are mutually 
dual [s^. These functions, plotted in Fig. 2, correspond to 

f \ ^ + y Inx-lny 2 

Cmin{x,y) ^ — , CKM[x,y) = , Cmax(a;,y) = — ■ — , (2.18) 

Zxy X — y X + y 

so the inverse, 1/c, is equal to the harmonic, logarithmic and arithmetic mean, respectively. It is known that among 
all normalised, self-transposed, operator monotone functions / there is a minimal and a maximal one [3^. They are 
distinguished above by the appropriate labels. The maximal concave function /max(^), the straight line t/2 + 1/2, 
induces the Bures metric and the Bures distance [s^l • 

If /(O) = 1/2 the metric induces the natural Riemannian metric on the subspace CP^^^ of pure states and is called 
Fubini-Study adjusted The Bures metric is the only monotone metric which is simultaneously Fisher-adjusted 

and Fubini-Study-adjusted. This property, valid for any matrix size A^, is easiest to discuss for N = 2. All metrics 
defined by p.l6|l are rotationally invariant, i.e. the uniform angular distribution (of eigenvectors of the density matrix) 
is independent of the radial distribution (of its eigenvalues). In the one-qubit case the metric depends only on the 
radius r = \/ + y^ + z'^ in the unit Bloch ball, and may be divided into the radial and the normal (tangential) 
components p^ . 



ds^ = -dr^ + g{ ]dn^ with g(t) = — - (2.19) 

The radial component is independent of the function /. It is easy to see that the value of f{t) aX t = {) determines 
the metric at the subspace of pure states, r = 1. If /(O) = 0, then the metric diverges at the manifold of pure states, 
which is the case for /km and /min- 

The relation between the Morozova-Chentsov function and the metric g = 1// shows that the larger the (normalized) 
function /, the smaller the metric and the corresponding distance. The Bures metric given by the maximal function 
/max(i), is therefore the minimal Fisher-adjusted, monotone Riemannian metric. These properties single out the 
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Bures metric and provide an additional motivation to study the geometry it induces in the set of mixed quantum 
states. 



III. BURES MEASURE 
A. Infinitesimal distance 



We are going to compute the Bures distance H2.4|l between two infinitesimally close density matrices of size N. Let 
us set PA ^ p, pb = P + Sp: 



pTpb pT =P + X + Y (3-1) 

where the matrix X is of order 1 in (Jp, while Y is of order 2. Squaring this equation we obtain to first and second 
order 

p^'^Sp p^/^ ^ Xp + pX , -X^^Yp + pY (3.2) 
or in the basis in which p is diagonal with eigenvalues > 

^1/2^1/2 ^ 

Xun = Sp^n ■ , Y^n = —{X )^M ■ . (3.3) 

Pi^ + Pfj. Pu + PtJ. 

Since Trp = 1, hence TiSp = and TrX ^ Tr Sp/2 0, while 

Thus we arrive at the resuh of Hiibner pjj for the Bures metric 

{ds)l = {D{p, p + 5p)f = i ■ (3-5) 

Note, that ii p^ = and p^ — 0, Spi^f^ does not appear and therefore terms where the denominator vanishes have to 
be excluded. 



B. The volume 

Let us rewrite the metric ()3.5|) distinguishing the diagonal elements 

u=l i'<^i + Pf" 

For diagonal matrices the second term vanishes and the above expression reduces to (I2.8|) . Since Trp = 1 so 
'^jj^i 5pyy = and not all Spui, are independent. Eliminating 5pNN we obtain 

V = V + ( y 5p..f— . (3.7) 

t^l p- t^l p- t^l p^ 

The metric (|3.7|l in the (X — l)~dimensional subspace is of the form gn^ = {5ik/ Pi + 1/ pn) with determinant Aetg — 
{pi+ P2 + + Pn)/{piP2---Pn) = 1/ det p using Trp = 1. Thus the volume element gains a factor y/detg — (detp)^^/^: 

1 2 
^Vb = T^W^idetp)-'/^ n n n dRc(p.^)^/2 dlm(p,^)) . (3.8) 
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The total volume of the set jv is given by: 



(3.9) 



with 6{p) = 6{pi)...6{pn) a step function that restricts integration to nonnegative p. We introduced the volume 
element 

N 

Dp^W dp,, W (V2 ARc{p,^)V2 dlm(p,^)) . (3.10) 

which differs from the usual definition by a factor 2^(^~i)/2. It follows from the metric for unrestricted Hermitian 
matrices: 

N 

(ds)2 = ^((5p..)2 + 2^|(5p.^|2 . (3.11) 

I/=l I/<jLi 

Above we considered the case of Hermitian density matrices (/? = 2). In the case of real density matrices (/3 = 1) 
there are no imaginary parts Ivup^n. 

C. Radial coordinates 

Let us diagonalize the matrices p 

p = UMJ-^ or p = OAC'-i (3.12) 
with U unitary for j3 = 2, O orthogonal for /? = 1, A = diag(pi, p2j Pn)- Then wc have 

dp = U[dK + U-^dUK- MJ-^dU]U-^ ov dp = 0{dh + 0-^AOK- K0-^dO\0-^ . (3.13) 
Let us first rewrite the Bures metric in a basis independent way, 

(ds)! = - / At Tr(dpe-''*dpe-''*) . (3.14) 
2 Jo 

This means it contains left and right multiplications with p. Inserting now the expressions for dp wc sec that again 
we can go to the representation where p is diagonal and obtain metrics in the cases of complex density matrices, 

(d^)^ = E ^ + E %-^l(t/-Mt/).,P , (3.15) 

and real density matrices 

(d^)' = E %^ + E %-^l(0-MO).,|^ . (3.16) 

In both cases the Bures measure has the following product form 

dVi'^) = dp(^\puP2,...,PN) X dvp. (3.17) 

The first factor depends only on the eigenvalues pi, while the latter, determining the distribution of the eigenvectors 
of the density matrix, is the unique, unitarily invariant, Haar measure on U{N) for /3 = 2 or on 0{N) for /? = 1. In 
the general case of complex density matrices, (/3 = 2), the Bures measure in the simplex of eigenvalues reads 

P.(A) = PsiPuP., ...,P.) = c j(P^ + P^ + - + P- z ll Yl iP^^, (3.18) 
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where the normalization constant Cn assures that the integral of P over all eigenvalues equals to unity. The measure 
was derived by Hall ^3 1 and was later analyzed in (3^ |39|| . To compute the area of the surface of the set of complex 
density matrices, and to find the volume of the set of real density matrices we will need a family of generalized 
constants Cn{oi,P) defined by 



1 



CN{a,l3) 



POO N , 



- Pn 



'I...N 

n 



{Pv - Pm)^ 



(3.19) 

Cn{1,2). They wiU be 



The standard Hall constants entering H3.18|l become a special case of this definition, Cn 
needed to compute the Bures volume of A^jv- 

Integrating the Bures measure over the eigenvectors distributed according to the Haar measure on U{N) or 0{N) 
we obtain 



d^]'') =2i-^[i+(^-i)^/2idetp-i/2 n n 



2(p. - P^.? 



Pu + Pfj. 



13/2 



(3.20) 



where vj^^^{N) denotes the volume of the flag manifold, Fl% = J7(7V)/[/(l)^ or Fl^ = 0(iV)/0(l)^ for /3 = 2, 1. 

■^N(N-i)/2 ^ (27r)^(^-i)/'^7r^/2 



J:' lag \- / 

Although these volumes 



and V^ll^iN) = Vol(F/« ) = 



(3.21) 



r(i)r(2)...r(7V) ^'^s' ' ^ r(i/2)r(2/2)..T(Af/2) 

were computed by Hua several years ago 40] , we provide their short derivation in To] and show that various results 
for the volumes of the unitary groups and flag manifolds present in the literature ,29, .41. .42„ 
choice of normalization. Formally we may write the volumes in the form: 



are due to different 



{13) 



N 



Flag 



{N) = \{{2^) 



(j-l)/3/2 



r(/3/2) 
r(j/3/2) 



(3.22) 



which shows that it is the product of the volumes of projective spaces with increasing dimension j ~ 1, real for /3 = 1, 
complex for /3 = 2, and the radius rescaled by a factor \f2 coming from the factor 2''^(^~i)/4 - see appendix of [l^ . 

Since the unitary or orthogonal rotations produce all permutations, in l|3.20|l we have to restrict to a certain order 
of eigenvalues, pi < p2 < ■■■ < Pn ■ In order to lift this condition we need to divide the result by the number of 
permutations, equal to A^!, 



(/3) 



2l-Ar[l+(Ar-l)/3/2] 



dpi...dpAr 

detpi/2 



<5(Trp-l) W 



(3.23) 



D. The surface 



The surface consists of the N pieces with one eigenvalue ~ 0. To calculate the metric on the surface (with pi = 0) 
we put first dpi = and then pi = (remember that terms with vanishing denominator have to be excluded): 



N 



4p, 



^9 2...N 

E 



(dp^)2 ^ 1^' (Pi. - Pm)^ 



N 



Pi' + Pt^ 



\{u-^duu\^ + Y,pMu-'au),, 



Then with the same reasoning as above the surface element is given by 



N-l 



dpz. 



j=2 Pv 



/2 



■2...^ 



N 



0/2 



and the total surface: 



qW _ q(/3) _ 9l-(Ar-l)(l+Ar/3/2) f TT dp,. 

J u=2Py 



dpi, 5{p2 + ... + Pat - 1) 



(A^-1)! 



^{pu - p,.y 



n 



Pu + Pt, 



<=2 



0/2 



(3.24) 



(3.25) 



<2^(7V), (3.26) 



where we have taken into account all permutations of the N — 1 non-zero eigenvalues. 
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E. The edges and hyperedges 

There are N{N — l)/2 edges which correspond to two vanishing eigenvalues. Hence there are {N — 2)1 different 
permutations of positive, generically different eigenvalues which form the edges. To calculate the metric on the edges 
we put first dpi = dp2 = and then pi = P2 = 0: 



N 



u=3 



3...N 

E 



Pu + P^i 



N 



\{U-'dUU\^ + J2 P,m-'dU)^,\' + \{U-'dUh,\') . (3.27) 



/i— 3 



We see that there do not appear terms {U ^dt/)^^ with v, p = 1, 2. Thus the corresponding coset space which has to 
be integrated over is the flag manifold U{N)/[U{2) x U{1)^-^] for /3 2 or 0{N)/[0{2) x 0(1)^-2] for /? = 1. Since 
their volume is equal to the ratio 'Vo\{Flf^)/Vol{Fl2), the edge element is then given by 



N-1 
i/=3 P'' 



dp_iy_ 

/2 



n 



<i(2) 



with the dimension dn 



(3.28) 



dn = {N ~n){l + {N + n- l)(i/2) - 1 



(3.29) 



The total area of the edges is 



c(/3) _ 9-<i2 



N 



v=3 Pv 



dpv 5(/03 + ■•■ + - 1) 



/2 



(A^-2)! 



.<.<M ^''^^^ K=3 



<L(2) 



(3.30) 



It is straightforward to generalize this to hyperedges: the manifold of states of rank N ~ n with n = 0, 1, ...A^ — 1 has 
the volume 



N 



n dpu S{p„+i + ... + piv - 1) 
-1/2 



2+1 Pf 



[N -ny. 



n+l...N ^ 

n n (2p« 



L'<IJ. 



Py + P^i 



K—n-\-l 



<l(«) 



(3.31) 



IV. BURES VOLUME 



Normalization constants 



The above formulae for the volumes of manifolds of states of rank N — n may be expressed in a compact way with 
use of the generalized Hall constants H3.19|l 



S 



N.n 



2-d„ + {N-n){N+n-l)l3/4 Vj^^^{N) 

{N~ny.CN-n{l + nP/2, (3) yp(fj^(„) 



(4.1) 



We are going to compute the constants Cat (a, (3) for integer a, (3 and then continue analytically to general positive a 
and (3. Let us go to variables pi = tf. Then the integral H3.19|l is an integral over the full sphere = tl+t2+...+t% — 1 
with an homogeneous function in R of order f3N{N — l)/2 + 2N{a — 1) as integrand. The radial integral may be 
replaced by a F function, which amounts to replacing S{R^ — 1) by e~^ . The result is the same up to a factor. 
Afterwards we may go back to the variables pi and obtain 



r[A^(2a + (TV - l)/?/2 - l)/2] 



N 



dpi, 



i/=i Pi' 



/2 



-(pi + ...+Pn) 



'1...N 

n 



jPu - Pp? 
Py + Pp 



0/2 



N 



(4.2) 



In the next step we want to eliminate the denominators in (|4.2(l . In this section we discuss the case /? = 2 of 
complex density matrices. 



1 



n 



(pi...piv)l/2 



Pm 



2—1 ^ i<j 



d KcXij d ImXij 



(4.3) 
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Using the volume element DX of a Hermitian Matrix X (|3.10|) we may rewrite this as 

TT ^ ^-^V2 2"^(A^-i)/2 [ DX e-Tr 



(4.4) 



This expression is basis independent and provides a way to describe the Bures probability distribution in a random 
matrix description. With this we may write for [3 = 2 



T(N^/2 + N{a-l)) ^ ^_^2 
Cjv(a,2) 



/2 2-N{N-l)/2_ 



Nl 



Now we use a formula which is an extension of the Ingham-Siegel integral ^ 

-N 



Dp e{p) det / DX e-T"- 



det 



S_ 
Tp 



6{p) = ANe{p) e 



-Tr tp 



(4.5) 



(4.6) 



valid for a positive Hermitian matrix e and with a constant Ajq. Let us derive this formula: 

det + Kp) = ^'""^ I d^,...d^N e-[^l(^/^''+^)'^i+-+<(^/^''+^)^~l 6{p) (4.7) 

with A^-dimensional complex vectors ipi, ...^ippf. The translation operator shifts the argument p of the (5-function by 
the matrix V'l ^ + ••• + ® i/'jv Then the integration can be done and yields 



det 



5p 



-N 



S{p) = e-T^ n-^ 



d^pi.. A^N S{p - (V'l <E)ipl + ... + ^/-AT (g) V'jv)) 



For p > we may rescale ip — > ^/pV'- The Jacobian cancels out and the result is equation H4.6|l . The constant An can 
be calculated from integrating this equation (|4.6|l over p. With partial integration we find 



2N{N~l)/2 



N 



An / Dp 0{p) e-^ 



det e 

This is valid for any Hermitian matrix e > 0. The coefficient is independent of the matrix e. To calculate An we 
put e = I. The integral can again be done using the Laguerre ensemble: 



1 

An 2^(^-i)/27V! 



N 



n r(i + j)r(j) = r(i)..T(7V) . 



(4.8) 



B. Main results 

Using equations (|4.51 14.61 14. 8() we may write 

r(7VV2 + Ar(a - 1)) _ 2-^(^-i)/2 iV! 



For a = 1 we obtain immediately by partial integration: 

r(ivV2) _ TT-^^/^m f DA 



Dp / DA detp° 



det ( — + 1 + 
Sp 



Cn{1,2) ANVj^l^iN) J det(l + A2) 



^-^V2 da 



A 



N 



Sip) . (4.9) 



(4.10) 



The last integral may be calculated in the complex plane and yields: (27r)^A^!2 ^ , such that the Bures-Hall constants 
are given by 



Cat — Cn{^, 2) — 



T{N^/2) 



7r^/2 r(i)...r(A^ + 1) 



(4.11) 
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This coincides with the known values for the Bures-Hall constants, C2 = 2/tt, C3 = SS/tt and C4 = and 
confirms some conjectures of Slater concerning C5 and Cg |3. The Bures volume of the set of complex density 
matrices of size N is now given by 

Ol-N(N+l)/2 1/(2) ( Af) 
. .(2) _ 2 ^ ^/ ^F\^.g[") _ „l_Ar2 n / 

^ " CAr(i,2) m ^ r(iV2/2) ^ ^ 

and for = 2 reduces to (|2.13|) . Note that the factor 2^~^ enters from the beginning in the definition. It disappears, 
if we increase ds (|3.6|l by a factor of 2. The Bures volume is exactly equal to the {N'^ — l)-dimensional volume of an 
{N'^ — 1) dimensional hemisphere with radius 1/2. Equations H4.11|l and H4.12|l constitute the central results of this 
work. 

C. Area of surface and edges 

For the surface and the edges of the Bures manifold we need a ^ 1. Scaling p — s- e"^/'^ p (T^l"^ for e any Hermitian 
matrix > shows that the integral H4.9|l can be obtained from 



Dpdet p° 



det I — + e 
op 



6 ^^-^ 



S{p) = Bn dete^-^-" (4.13) 



with Bjs[ independent of e. To compute B]\j we put e = I and go back using equation (|4.6|) . This again leads to an 
integral over the Laguerre ensemble and we obtain: 

T/(2) (nr^ N ^ -Tf ^ ■ -t\ 

j=l J=l 

Such an integral appeared already for the H-S metric 19]. Thus we have 

r[iVV2 + 7V(a - 1)] 7r-^'/2iV! f DX 



C^(a,2) AmV^I^{N) J det(l+X2)Af+"-i 



r(i/2) + 



Thus the ensemble is reduced to a kind of the Lorentz ensemble. The last integral can be found in the book by Mehta 
[i^ . (page 350) following from the Selberg's integral: 

J(l, 1, iV + a - 1, iV + a - 1, 1, iV) = ^N^-N^+N-2Nic^i) TT r(l+j)r(Ar + 2a-- l^j) 

fJl r(iv + a- iY 

As a result we obtain 

T{N^/2 + N{a^l)) ^ «/2„_^^+^_2^(„_i) r(l+j)r(2a^2 + j) 

CN{a,2) n r(a + ,-l) ■ ^^-^'^ 

From equation H4.1|l we obtain for the volumes of submanifolds of rank — n of the space of complex density matrices 

2-.„ + (Ar-n)(Ar+»-l)/2 y^f4(iV) 

CA._„(l + n,2)(iV-n)! y^^i^^^ ' ^^''"^ 

The result is rather simple [s^ 



r[(d„ + l)/2] 11 [(j+n)!]2 ^^-^'^ 
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and reduces to H4.12|l for n = 0. Here (i„ = N"^ — — 1 denotes the dimensionality of each submanifold. Observe that 
the above prefactor equals to the volume of a (i„-dimensional hyper- hemisphere with radius 1/2. The above results 
allows us to find the ratio between the surface and the volume, 

^2)^ 0Fr(iVV2-l/2)''- ^^-'"^ 

It behaves as N'^ and in the limit of large matrix size it grows proportionally to the dimensionality N"^ — 1 of the 
space of mixed states. Note that such a ratio for the hemispheres of dimensionality N"^ is proportional to N . For 
pure states n = — 1 we obtain the volume 

S^S]n-i = Vol(CP^-i) = ^^^-^ , (4.21) 

which is equal to the volume of a 2{N — 1) dimensional real ball of radius 1 and this equals the volume of a (iV — 1)- 
dimensional complex projective space, which is a correct consequence of the Bures metric for pure states. By the way, 
up to a scale factor this is the same for the Hilbert-Schmidt metric 19]. Differences come only for mixed states. 

V. REAL DENSITY MATRICES & LORENTZ TYPE ENSEMBLE 

To compute the volume of the set of real density matrices we need to find the generalized Hall constants (|3.19|) for 
(3—1. To do so we write 

(pi...pw)i/2 11 (p,+p^)i/2 J 11^11^ 
Let us introduce the volume element of a real symmetric matrix X with the metric 

N 

{dsf = Tr(dX)2 = (d^^O' + 2 '^^l ' (5-2) 
DX^\[dXu \{{dXjkV^). (5.3) 

i j<k 



Then we have 



With this we may write the formula for /3 — I 



i...pjv)i/2 ,J^Mp.+p^)i/2 J 



(5.4) 



r(jV(a-l/2 +(jy-l)/2)) ^ (2^)-Ar(^+i)/4 r^p^,„Ap^ rr \p^- p^\Ueip'^-' f DAb-^^Mi+a^^ . (5.5) 



Now with e a real symmetric positive matrix we want to calculate for real symmetric p\ 

-N/2 



det 



= e-^''pTT-^" f d4>i...d^N 6ip-i4>i®^Pj + ... + ^n®^n)) , (5.6) 



with real vectors V'l, V'JV- For p > we may rescale ip — > y/pi^, then the volume element is cx det p^/^ and the 
(5-function oc det This implies 



det ( ^ + e 

Sp 



-N/2 

Sip) = Kn e{p) e-^ det (5.7) 
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with Kj^i a constant independent of e. With partial integration we find 

r)JV(Ar-l)/4 r 

Kn Dp 6{p) e-^' detp-i/" 



det 

To calculate Kn we put e = I and find again with the help of the Laguerre ensemble 01 (page 354) 



1 VF^/ig(^^) 

Kn 2^(^-i)/4iV! 



N 

n 



r(l+j/3/2)r[l/2+(j-l)/3/2] 

r(i + /3/2) 



(5.8) 



To obtain the generalized normalization constants for (3 = 1 and arbitrary a we may write 



r[A^(2a-l + (A^-l)/2)/2] _ (^2TT)-^^^~^y^2^/^Nl 



Cjv(a,l) 



Dp / DX detp"-i/2 



det I — + 1 + 



Scaling p ^ e p e for e any real symmetric matrix > shows the identity 



L»pdet p° 



det I — + e 

op 



-N/2 



5{p) = Ln dete(i-^-2a)/2 



5{p) . 
(5.9) 

(5.10) 



with Ln & constant independent of the matrix e. To compute Ln we put e = I and go back. Again we obtain an 
integral over the Laguerre ensemble with the result: 



r,JV(W-i)/4 TT r[a + (j - l)/2] 



/2 + (j-l)/2] 



(5.11) 



This formula is very similar to (|4.14() for the constant Bn in the complex case. In fact it is not difficult to rewrite 
both results in one expression as a function of 13. This allows us to write down an equation, which is valid as well in 
the complex or in the real case, for (3 = 2 and (3=1, respectively, 



T[N{2a -1 + {N- l)/3/2)/2] ^ T[a + [j - l)/3/2] 



CN{a,P) 



r(i/2) 



+ 00 



N 



n n+.2)(^-iw2+^ n - -^1' ■ (5-12) 



As result we have a Lorentz type ensemble. The last integral can be found in the book by Mehta (page 350) with 
7 = {N ~l)(3/2 + a: 



2^) \i r[l + U + l)(3/2] r[27 -{N + j- l)/3/2 - 1] 
r(l + /3/2) [r(7-j/3/2)]2 



J(l, 1, 7, 7, P/2, N) = -^N^NH,N-l)p/2+2~2^) -Q 



(5.13) 



Inserting this into eau. H5.12|l we obtain 

r[7V(2Q-l + (iV-l)/3/2)/2] 
CN{a,(3) 



rN/2 



N 



2Af((JV-l)/3/2+2(Q-l)) 



r(l+j/3/2) r[(JV - j)/3/2 + 2a - 1] 
r(l + /3/2) T[{N~3)p/2 + a] 



(5.14) 



From the above equality expression we may read out the explicit form of the generalized Hall constants Cjv(a,/3). 
They allow us to obtain a fairly general expression for the Bures volume of the submanifold of the states of rank 
iV — n of the set of complex {(3 = 2) or real {[3=1) density matrices 



S 



N,n 



2_,„ r(j/3/2) r[l + (2n + j-l)/?/2] 

r((d„ + l)/2) f \ r[(n + j)P/2] r[l + (n + j - l)(3/2] 



(5.15) 



where d„ = (TV — n)[l + {N + n ~ l)(3/2] — 1 represents the dimensionality of the manifold. In a sense, this formula 
may be considered as a brief summary of the paper: for n = the last factor simply equals unity and (|5.15l) gives the 
Bures volume of the entire space of density matrices, equal to that of a do-dimensional hyper-hemisphere with radius 
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1/2. In the case n = 1 we obtain the volume of the surface of this set, while for n = iV — 1 we get the volume of the 
set of pure states 

^N,N-i- r(iV/3/2) ' ^^-^^^ 

which for P = 1(2) gives correctly the volume of the real (complex) projective space of dimension A'^ — 1. This result 
can be seen directly from equations (|3.31I3.22() . In the case of a real projective space the volume is equal to that of a 
hemisphere of dimension — 1 and radius 1: tt^^^ /r{N/2). The final formulas are probably correctly continued to 
general positive though we have strictly proven them only for /3 = 1, 2. 



VI. CONCLUDING REMARKS 



In this work we have computed the volume of the set of complex (real) density matrices with respect to the Bures 
measure. Interestingly, this volume is equal to that of a hyper-hemisphere of radius 1/2 and the same dimensionality. 
However, this fact does not mean that the geometry of the set of complex density matrices induced by the Bures 
metric is identical with the geometry of a hemisphere of dimensionality iV^ — 1 embedded in M.^ . This is the case 
for iV = 2, for which the set A^2 is isomorphic with the Uhlmann hemisphere ^S'^, of a constant constant positive 
curvature, embedded in E^. 

On the other hand, in the general case, N > 3, the curvature of the set A^2 induced by the Bures metric is not 
constant. To explain this fact, rigorously analyzed by Dittmann ,42j, consider the set of pure states for = 3. This 
manifold forms a complex projective space, CP^, and contains, on one hand, the sphere — CP^ and on the other, 
the real projective space KP^. The latter is equivalent to a hemisphere, (with all antipodal points identified), so it 
has a different geodesic length and curvature different than CP^ . Hence the set of mixed states for N > 2 analyzed 
with respect to the Bures metric, which preserves the geometry of the manifold of pure states, cannot have a constant 
curvature. Note that the Hilbert-Schmidt distance induces a flat geometry of the space of mixed states, (i.e. with 
the curvature equal to zero), for an arbitrary dimensionality N. 

The Bures measure, related to statistical distance. Fisher information and Jeffreys' prior, is often used to define and 
compute the probability that a random density matrix satisfies a certain condition. For instance such an approach was 
applied to find the probability that a random state acting on a given composed Hilbert space is separable [43, Ha, |43 . 
In this work we computed the Bures-Hall normalization constants, which may be useful to extend such calculations 
for higher-dimensional problems. 

Results presented in this work may also be considered as a contribution to the theory of random matrices. A closest 
work in this field is due to Eynard and Kristjansen who analyzed the 0(n) model on a random lattice. We 
introduced a different class of ensembles of random matrices defined by the joint probability density of eigenvalues 

P^'''{Pi,P2,-,Pn) = Cn[1,I3)— , II 7 ■ (6.1 

with normalization constants computed in H5.14|l . which for (3 = 2 and (3=1 gives the Bures measure on the set 
of complex and real density matrices, respectively. The density of eigenvalues for these ensembles will be studied in 
a forthcoming publication [sj- It might be also instructive to analyze a family of such ensembles by allowing the 
parameter (3 to admit arbitrary real values. In the limit /3 — > this distribution tends to the Dirichlet distribution of 
order 1/2. It describes the distribution of squared components of a real random vector, i.e. a column (or a row) of a 
random ortogonal matrix generated according to the Haar measure on 0{N) fs^fs^l. 
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